Can Ballistic Electrons Probe the Electronic Spectra of Individual Buried Molecules? 
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A theoretical study is presented of the baUistic electron emission spectra (BEES) of individual 
insulating and conducting organic molecules chemisorbed on a silicon substrate and buried under 
a thin gold film. It is predicted that ballistic electrons injected into the gold film from a scanning 
tunneling microscope tip should be transmitted so weakly to the silicon substrate by alkane molecules 
of moderate length (decane, hexane) and their thiolates that individual buried molecules of this 
type will be difficult to detect in BEES experiments. However, resonant transmission by molecules 
containing unsaturated C-C bonds or aromatic rings is predicted to be strong enough for BEES 
spectra of individual buried molecules of these types types to be measured. Calculated BEES 
spectra of molecules of both types are presented and the effects of some simple interstitial and 
substitutional gold defects that may occur in molecular films are also briefly discussed. 

PACS numbers: 81.07.Nb, 73.63.-b, 81.07.Pr, 73. 23. Ad 



I. INTRODUCTION 

During the last decade there has been growing interest 
in molecular electronics, stimulated largely by the ex- 
perimental realization of molecular wiresi^i^, systems in 
which a single organic molecule or a few molecules carry 
an electric current between a pair of metal contacts. Hy- 
brid molecule/silicon nanoelectronic devices are another 
intriguing possibility and research that may lead to their 
creation is also being pursuedi^i^i^iii^i^ii^iiiii^ For many 
potential nanoelectronic applications it will be neces- 
sary to sandwich a molecular layer between metal and/or 
semiconductor layers that act as contacts and also protect 
the molecular layer from the environment. However, this 
makes the molecules inaccessible to direct nanoscopic ex- 
perimental probes such as the scanning tunneling micro- 
scope (STM). Thus the structures and electronic proper- 
ties of such buried molecular layers are, to a large degree, 
still the subject of conjecture. 

Ballistic electron emission microscopy (BEEM) is a 
technique that can probe the nanoscale spatial structure 
and electronic spectra of buried interfaces F^^i^^ In BEEM 
experiments a bias voltage is applied between an STM tip 
and a thin metal overlayer covering a buried substrate, 
usually a semiconductor, which is kept at almost the 
same electrochemical potential as the metal overlayer^iS 
Electrons are injected from the STM tip into the metal 
layer and travel ballistically through the metal. Some of 
them are transmitted through the metal/semiconductor 
interface into the semiconductor substrate. The electric 
current transmitted through the interface is measured 
as a function of both the lateral position of the STM 
tip and the bias voltage applied between the tip and 
metal overlayer. This allows the imaging of defect struc- 
tures at the interface and spatially-resolved spectroscopy 
of their electronic structure, with nanometer resolution. 
Recently BEEM experiments have been initiated on sys- 
tems in which self-assembled molecular monolayers are 
present between the metal overlayer and semiconductor, 
specifically, alkane dithiolates between GaAs and gold^i^ 



Since in such BEEM experiments a significant bias is not 
applied across the molecular layer, BEEM spectroscopy 
(BEES) of buried molecules (unlike the usual spectro- 
scopies based on measurements of current-voltage charac- 
teristics of molecular diodeaii^i^) is not subject to compli- 
cations due to bias-induced charging effects and electric 
fields in or near the molecular layer. Thus in addition to 
being a unique probe of buried structures with nanometer 
spatial resolution, BEEM can provide important spectro- 
scopic information relevant to molecular electronics that 
is complementary to that obtained from measurements 
of diode current-voltage characteristics. However, be- 
cause of the absence to date of any appropriate theory, 
it has been difficult to interpret the results of BEEM 
experiments for metal/molecule/semiconductor systems 
unambiguously: For example, in the recent experimen- 
tal study^^ a BEEM signal could only be detected in 
a few small, isolated patches of the sample and it has 
been unclear whether this should be interpreted as evi- 
dence that most of the sample was occupied by an insu- 
lating alkane-dithiolate self-assembled monolayer (SAM) 
that did not transmit a detectable BEEM signal and that 
the observed signal was coming from occasional patches 
of defects in this SAM, or that to the contrary the de- 
tected BEEM signal was being transmitted through small 
patches of alkane-dithiolate SAM surrounded by much 
larger regions of a more strongly insulating material (such 
as alkane-dithiol multilayers) between the semiconductor 
and the gold overlayer. 

In this article a theoretical study is presented of ballis- 
tic electron emission microscopy/spectroscopy of single 
organic molecules between a thin gold overlayer and a 
silicon substrate. It is predicted that alkane molecules of 
moderate size (hexane and decane are studied) attached 
to the Si substrate by a C-Si bond (and whether thiol- 
bonded or not chemically bonded to the gold overlayer) 
should transmit the ballistic electrons so weakly that such 
single molecules are expected to be below the threshold 
of detection for conventional BEEM equipment. How- 
ever, related single molecules containing unsaturated C- 
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C bonds or aromatic rings are predicted to exhibit suffi- 
ciently strong transmission resonances that their obser- 
vation in BEEM experiments should be feasible. Thus it 
is predicted that single molecules with unsaturated C-C 
bonds or aromatic rings embedded in an alkane or alkane 
thiolate SAM between a silicon substrate and gold over- 
layer should be detectable in BEEM experiments with 
present day equipment. The calculated BEEM spectra 
of the molecules are also presented and the effects of 
some simple interstitial and substitutional gold defects 
that may occur in molecular films are also briefly dis- 
cussed. 

This paper is organized as follows: The theoretical 
model of the structures and electronic structures of the 
systems studied and the approach used to calculate the 
BEEM currents in these systems are described in Section 
im The results obtained are presented together with their 
interpretation and implications for experiment in Section 
mil The main conclusions are summarized in Section ITVl 



II. MODEL AND THEORETICAL APPROACH 

A. Extended Molecule 

As in other transport calculations in molecular 
electronics^iiSii^iiiii^ii^i^ai^iS^ the theoretical model 
studied is based on an extended molecule that includes 
the molecule itself as well as clusters of nearby atoms be- 
longing to the electrodes that transmit electrons to and 
from the molecule. In this way the chemical bonding 
between the molecule and the electrodes is taken into ac- 
count as well as the electronic structures of the molecule 
and the electrodes. In the present systems there are three 
electrodes: the STM tip, the metal layer between the tip 
and the molecule, and the silicon substrate. An atomic 
cluster from each of these is included in the extended 
molecule, as is shown for an ethylmethylbenzene molecule 
in Fig[l] A cluster of 10 Au atoms represents the (111)- 
oriented monoatomic Au tip. The Au film between the 
tip and molecule is represented by a 591 Au atom clus- 
ter that is approximately cylindrical and presents (111)- 
oriented facets to the STM tip and the molecule. A hemi- 
spherical 390 Si atom cluster with the danghng Si bonds 
passivated by H atoms and its flat (111) surface facing 
the molecule represents the Si substrate. The molecules 
studied here bond to the Si substrate through a single 
covalent C-Si bond. For those molecules that are thiol- 
bonded to the gold film, the sulfur atom is positioned 
'--^ 2. 2 A below a hollow site between three gold surface 
atoms. Density functional theory was used to estimate 
the atomic geometries of the molecules and of their bond- 
ing to the semiconductor substrate and metal overlayer.^^ 
However in cases where the end of the molecule adjacent 
to the gold layer terminates in a CH3 group (which does 
not bond to the gold), the distance between the gold 
surface plane and the closest H atom of the molecule is 
somewhat arbitrary, and was assumed to be SA. 



B. Electronic Structure 

Most theoretical work on electronic transport in 
molecular wires with metal contacts has been based 
on semi-empirical tight-binding models or Kohn- 
Sham density functional calculations of the electronic 
structure i'^'-^''i-^^'-^^i^°'^^i^^ It was found in the present 
study that, in order to realistically model ballistic trans- 
mission of electrons through the gold film between the 
STM tip and molecule in BEEM experiments, several 
hundred Au atoms need to be included in the cluster 
that represents the gold film in extended molecule. This 
makes the density functional approach impractical for 
the transport problem. Furthermore previous theoret- 
ical work modeling the current-voltage characteristics 
of molecules on silicon found Kohn-Sham density func- 
tional calculations to be unsuitable for treating the elec- 
tronic structure of the silicon substrate, and to yield re- 
sults that disagree qualitatively with experimental STM 
datafi^ However, transport calculations based on semi- 
empirical tight-binding models can be performed for sys- 
tems with large numbers of atoms and have successfully 
explained experimental current-voltage characteristics of 
molecules connecting gold electrodea^ i^^'^^ and of molec- 
ular wires on silicon substratesJ^ Thus a hybrid tight- 
binding model of the electronic structure that combines 
parameters derived from semi-empirical considerations of 
quantum chemistry and from ab initio band structure cal- 
culations is adopted here. As is explained below, a num- 
ber of parameters that enter the model are not known 
accurately, and in order to proceed further it is neces- 
sary to make plausible but not rigorous approximations 
and assumptions regarding their values. The sensitivity 
of the present findings to the values of these parameters 
was explored in the course of this work and it was found 
that the qualitative conclusions reported here are robust 
to reasonable variations in parameter values. The reader 
is referred to previous paper a^^'^'^i^^ for extended discus- 
sions of the reliability of such techniques. 

The semi-empirical extended Hiickel model of quantum 
chemistrySIiSi, is used to describe the electronic structure 
of the molecule and coupling between the molecule and 
silicon substrate. However in order to model the elec- 
tronic structure of the silicon substrate satisfactorily, the 
extended Hiickel model requires modification and this is 
done here as in Refill: Extended Hiickel theory describes 
molecular systems in terms of a small set of Slater-type 
atomic orbitals {|(/>i)}, their overlaps Sij = {(/)i\(l>j) and 
a Hamiltonian matrix Hij = The diagonal 

Hamiltonian elements Ha = are the atomic orbital 
ionization energies. Non-diagonal elements are approxi- 
mated by 

H,, ^ KS,,ie, + e,)/2 (1) 

where K is a, phenomenological parameter usually cho- 
sen to be 1.75 for consistency with experimental molec- 
ular electronic structure data. However, in the present 
work, as is discussed in detail in Refll2l, for Si atoms 
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the single parameter K is replaced by a set of param- 
eters Kij whose values are chosen so that the modified 
extended Hiickel model obtained in this way yields an 
accurate description of the band structure of crystalline 
silicon. The electronic structures of the gold film and 
STM tip are also described by a tight-binding model with 
a non-orthogonal basis but in this case the tight-binding 
parameters Hij and Sij from Ref. [2^ that are based on 
fits to ab initio calculations of the electronic structure of 
gold are used. The atomic orbital energies Ha of Ref. 
|29| for gold are defined up to an arbitrary additive con- 
stant which is chosen in the present work so as to yield 
reasonable values of the energy offsets Am.Au between 
the HOMO levels of the molecules M studied and the 
Fermi level of the gold film. The values of the offsets 
Ajv/jAti are not known accurately at present either theo- 
retically or experimentally, their determination being an 
important open problem of molecular electronics. For the 
present purpose they are estimated as the differences be- 
tween the work function of Au and the HOMO energies 
of the respective isolated molecules obtained from den- 
sity functional theory, an approximation that has been 
used successfully to align the Fermi energy of gold with 
the HOMO of benzenedithiolate^^ and has also been ap- 
plied to a variety of other metal-molecule junctions.— 
The effects of the potential difference between the STM 
tip and the gold film on the electronic structure of the 
system are included in the present model by offsetting 
the atomic orbital energies Ha for gold STM tip up- 
wards in energy from those of the gold film by an amount 
\eV\ where V is bias voltage between the tip and Au 
film. The value of the energy offset /^si,Au of the sili- 
con conduction band edge (at the surface of the silicon 
substrate) above the Fermi level of the gold film is an- 
other important parameter that is not known accurately 
for gold/molecule/silicon heterojunctions. Here it is as- 
sumed, for simplicity, that l^si,Au = 0.84ey, a typical 
experimental value of the Schottky barrier measured in 
BEEM experiments on Au/Si(lll) interfaces)^ the Si 
orbital energies Ha = ti are adjusted (all equally) so 
as to place the LUMO level of the Si cluster of the ex- 
tended molecule 0.84ey above the Fermi level of the gold 
film. Finally, the Hamiltonian matrix elements between 
the STM tip and gold film and between the gold film 
and molecule are estimated from Eq.([T|) with K = 1.75. 
Because of the non-orthogonality of the basis states em- 
ployed in the present model, the above adjustments of the 
diagonal matrix elements Ha of the tight binding Hamil- 
tonian require-^- corresponding adjustments AHij of 
the non-diagonal elements Hij , which were taken to be 



AH,, = S^j{E,+Ej)/2 (2) 



where Ek is the shift of the diagonal element Hkk- 



C. Transport 

In a BEEM experiment the electron flux that passes 
from the STM tip to the gold film subsequently separates: 
A small part of it (that constitutes the BEEM current 
^beem) continues through the molecule to the Si sub- 
strate and on to the electron drain electrode, while most 
of it (referred to henceforth as the film current /film) 
proceeds laterally through the gold film directly to the 
drain electrode. The total current /tip passing through 
the STM tip is thus given by /tip = /beem + /film- 

The present calculations of /beem and /film are based 
on Landauer-Biittiker theory.— For BEEM experiments 
the gold film, the Si substrate and the electron drain 
can all be considered to be at the same electrochemical 
potentialji^ Also, in practice, electron transmission from 
the source to drain via the gold film, molecule and silicon 
substrate, and directly via the gold film are mutually in- 
coherent processes. Therefore it follows from Landauer- 
Biittiker theory that the currents /beem and /film at a 
given value of the applied source-drain bias voltage V are 
related to the multichannel electron transmission proba- 
bilities Tx{E, V) at energy E from the source electrode to 
the drain via x where x can be the BEEM path through 
the gold film, molecule and Si substrate (x — BEEM), 
or the direct path through the gold film to the drain 
{x = FILM), according to 

LAV)^~- j dET,{E,V){f{E,fi,)-f{E,fia)) 

(3) 

where f{E, fii) = l/{exp[{E — ^i)/kT] + 1) and is the 
electrochemical potential of the source (z = s) or drain 
(i — d) electrode;^? 

The source and drain electrodes are modelled here as 
arrays of ideal single-channel leads: One such drain lead 
is coupled to the Is orbital of each of the hydrogen atoms 
that passivate the dangling bonds terminating the Si clus- 
ter of the extended molecule, and one is coupled to each 
p valence orbital of the surface atoms of the Si cluster 
(except for those of the H and Si atoms in the vicinity of 
the Si atom to which the molecule bonds), a total of 675 
ideal leads. An ideal drain lead is also coupled to each 
s, p and d valence orbital of the surface atoms of the Au 
cluster representing the the Au film, except for those Au 
atoms in the immediate vicinity of the molecule and STM 
tip, totalling 2709 leads. Similarly an ideal source lead is 
coupled to each valence s, p and d orbital of the surface 
atoms of the cluster representing the STM tip, except for 
the terminal atom of the tip, a total of 81 leads. Each 
ideal lead k was modelled as a semi-infinite tight-binding 
chain with a single orbital per site, a site energy ak and 
nearest neighbor hopping matrix element (3. ak was cho- 
sen to be equal to the energy (including the shift, if any, 
due to the applied bias) of the hydrogen, silicon or gold 
orbital to which lead k was coupled. The hopping matrix 
elements (3 were all taken to have a magnitude of 6 eV, 
sufficiently large that the eigenmodes of all of the ideal 
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leads in the energy window between fig and fid in eq. ([3]) 
be composed of counter-propagating (left and right mov- 
ing) waves. The coupling matrix elements Wk between 
the ideal leads and their respective hydrogen, silicon and 
gold orbitals were also set to (3. 

As well as mimicking macroscopic electrodes by trans- 
mitting an ample electron flux to and from the sys- 
tem, the above large numbers of ideal source and drain 
leads have a similar effect to phase-randomizing Biittiker 
probes^ in minimizing the influence of dimensional res- 
onances due to the finite sizes of the gold and silicon 
clusters employed in the model; it was verified that the 
ideal leads described above meet these requirements by 
comparing the results obtained in test calculations with 
differing numbers of ideal leads, differing values of the 
lead parameters, and differing numbers of atoms in the 
Au and Si clusters of the extended molecules. Calcula- 
tions of the currents /beem and /film for cylindrical Au 
atom clusters with the same height (1.4nm) as that of the 
591 Au atom cluster in FigU] but with differing diame- 
ters confirmed that the diameter of this cluster (approx- 
imately 2.5nm) was sufficient for these currents to have 
converged with increasing cluster diameter, i.e., that elec- 
tron transport through this finite cluster adequately rep- 
resents ballistic propagation of electrons through a thin 
but laterally extended gold film. 

To calculate Tx{E,V) and hence evaluate eq.Q, the 
transformation to the alternate Hilbert space described 
in Refs. 'sFand^^ was made, mapping the non-orthogonal 
basis of atomic orbitals to an orthogonal basis. The 
Lippmann-Schwinger equation 



|*fc) = |$o,fe)+Go(/^)M^|*fc) 



(4) 



describing electron scattering between the source and 
drain leads was solved numerically for j^'fc) in the al- 
ternate Hilbert space.— In eq.(j4]) Go{E) is the Green's 
function for the decoupled system (i.e., with the cou- 
pling between the ideal leads and the extended molecule 
switched off), |$o,fc) is the eigenstate of the decoupled 
ideal source lead k with energy E and is the cor- 
responding scattering eigenstate of the complete system 
with the coupling W between the ideal leads and the 
extended molecule switched on. The scattering ampli- 
tudes tjk from the ideal source lead k to drain lead j at 
energy E were extracted from the scattering eigenstates 
\^k) and the transmission probabilities that enter eq. ([3|) 
were then calculated from 



\tik\ 



(5) 



where Vk and Vj are the electron group velocities in ideal 
leads k and j respectively at energy E. The sum over j in 
Eq. ([5]) is over those drain leads connected to electrode 
X, i.e., the gold metal film or the Si substrate as discussed 
above, while that over k is over the leads connected to 
the STM tip. 



III. RESULTS 

For the terminal atom of the gold STM tip over an 
atom at the center of the top surface of the Au cluster 
that represents the gold film in FigUl and separated from 
that surface atom by a nearest neighbor distance of gold 
(2.88A), the STM tip and gold film are connected by an 
atomic quantum point contact. For this geometry the 
present calculations of the electron transmission proba- 
bility T'film(-E', 0) between the tip and gold film at zero 
applied bias yield values close to 1 (per spin) for ener- 
gies E around the Fermi level of gold, consistent with 
previous experimental^i^^ and theoretical^^s^S studies of 
gold atomic point contacts. However, in BEEM exper- 
iments the STM tip is normally further from the gold 
film, at distances in the STM tunneling regime. Un- 
less stated otherwise, the results of the BEEM calcu- 
lations to be presented here will be for such a larger 
separation (4.33A) between the tip atom and the gold 
film for which the calculated the electron transmission 
probability Tfilm(/'_f, 0) between the tip and gold film 
at the gold Fermi energy and zero applied bias is ap- 
proximately 10~^. The calculated electron transmission 
probabilities T-q-e-em{E ,V) from the electron source to 
the drain via the STM tip, gold film, molecule and sili- 
con substrate are shown in Fig. [5] for two representative 
molecules that are shown in the insets. The solid (dot- 
ted) curves are for zero (2V) bias V between the STM 
tip and gold film. The Fermi level of the gold film is at 
E=0 in all cases, and the silicon conduction band edge at 
the silicon surface Ec is indicated by the vertical dashed 
line. The lower plot is for decanethiolate (C10H20S) thiol 
bonded to the gold film below a hollow site between 3 gold 
atoms and attached to the silicon substrate by a single 
carbon-silicon bond. This molecule is insulating, having 
a large HOMO-LUMO gap, and thus its transmission in 
a BEEM experiment is predicted to be very small result- 
ing in T'beem(-£', ^) < I0~^^ in the range shown. The 
upper plot is for another chain molecule (CioHigS) also 
with 10 carbon atoms and bonded similarly to the gold 
film and Si substrate, but in this case the 5th and 6th 
carbon atoms each have one H atom attached instead of 
two so that there is a double bond connecting these two C 
atoms; see upper inset. This double C-C bond gives rise 
to a molecular state near E = 1.65eV where the trans- 
mission T'beem(-E', y) is greatly enhanced (at its peak by 
almost seven orders of magnitude relative to that for the 
decanethiolate) although it is still very weak (< 10^^) in 
absolute terms, as expected for a transmitting state asso- 
ciated with a double bond surrounded by large potential 
barriers. However, there is no other resonantly transmit- 
ting state in this system in the energy range above and 
reasonably close to the Fermi level of the gold film, the 
range relevant to BEEM experiments. Therefore despite 
its weakness this resonance is the dominant BEEM trans- 
port mechanism for this system as will be seen below. 
Notice that the application of a 2V bias does not change 
the transmission probabilities Tbeem(-£', V) by much on 
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the scale of Fig. [2] (less than a factor of 2) because the 
main transport bottleneck here is the molecule and the 
application of a bias between the tip and gold film does 
not affect the molecule significantly. 

In BEEM experiments the distance between the tip 
and gold film is not known accurately, as is the case for 
most STM experiments where only changes in the tip- 
surface separation are measured. However, the transmis- 
sion probabilities Tbeem and Tfilm and currents /beem 
and /film depend strongly on the tip-surface separa- 
tion. This makes a direct comparison between calculated 
and experimental values of these quantities problematic. 
However, it is intuitively plausible that /beem and /film 
should scale in the same way with the tip-surface sepa- 
ration in the tunneling regime and thus /beem //film 
should be nearly independent of the tip-surface separa- 
tion. This has been confirmed by calculations carried 
out in the present study. Thus the ratio /beem //film 
lends itself better than /beem or /film to comparison 
between theory and experiment and it will therefore be 
the focus of attention in what follows. Another advan- 
tage of studying this quantity theoretically is that spec- 
troscopic BEEM (i.e. BEES) measurements of /beem 
vs. the bias voltage applied between the STM tip and 
metal film are normally taken not at constant tip-surface 
separation but at constant tip current /tip, which, for 
the systems studied here, differs negligibly from the film 
current /film- Therefore because /beem //film is ap- 
proximately independent of the tip-surface separation, 
plots of /beem//film vs. bias voltage should differ from 
experimental plots of /beem vs. bias by only a scale fac- 
tor that is independent of the bias, making qualitative 
comparisons between calculations and BEES current- 
voltage characteristics appropriate. A typical experimen- 
tal BEEM setup can readily detect BEEM currents with 
current ratios /BEEM//tip greater than roughly 10~^; this 
number will be referred to below as the "nominal BEEM 
sensitivity threshold" (NBT). 

In the present calculations the silicon cluster that rep- 
resents the silicon substrate in the extended molecule is 
of necessity small even though it includes hundreds of 
atoms; the silicon hemisphere shown in Fig [1] has a radius 
slightly less than 1.5nm. Because of its small size it trans- 
mits electrons appreciably by quantum tunneling even at 
energies below the conduction band edge Ec in Fig [2] 
where there are no eigenstates of the extended molecule 
with a strong silicon content. Because of this the cal- 
culated transmission Tbeem (/i', V) shown in in Fig[21 al- 
though it is extremely weak below Ec, does not decline 
all the way to zero for E < Ec- Thus while the present 
model should describe the BEEM electron transmission 
reasonably well at energies higher than Ec, it does not 
describe the details of the onset of BEEM current at the 
threshold voltage V — Ec accurately for macroscopic 
silicon substrates. However, as will be seen below, for 
the molecules studied here, the current ratios /BEEM//tip 
in the near-threshold regime are very far below nominal 
BEEM sensitivity threshold, and detailed study of this 



regime is therefore left for future work. Here, for sim- 
plicity, in calculating /beem from Eq. Tbeem{E, V) 
will be set to zero for E < Ec- 

The calculated current ratios /BEEM//tip for some rep- 
resentative saturated chain molecules are shown in Fig. 
[31 C10H21 is decane bonded to the silicon substrate 
through a single C-Si bond and not bonded chemically 
to the gold film. Even at an STM tip bias of 4V the cal- 
culated /beem //tip for this molecule is 6 orders of mag- 
nitude below the nominal BEEM sensitivity threshold. 
For C10H20S (decanethiolate depicted in the lower inset 
of Fig. [2]) the thiol bond to the Au film results in pre- 
dicted BEEM currents two orders of magnitude larger 
than for C10H21. However, the calculated current ratios 
/BEEM//tip for both C10H21 and C10H20S are still so low 
as to be well below the NBT not only for single molecules 
but even for self-assembled molecular monolayers where 
hundreds of molecules may be contributing together to 
the BEEM current. For C6H12S (hexanethiolate which 
is similar to decanethiolate but with fewer carbon atoms 
and so presents a thinner tunnel barrier than decanethio- 
late) the calculated current ratios /beem //tip are two or- 
ders of magnitude higher than for decanethiolate but still 
well below the NBT. Thus detecting a BEEM signal from 
single alkane or alkanethiolate molecules on silicon or 
even from self-assembled monolayers of these molecules 
is expected to be challenging experimentally except for 
quite short molecules of these types. This is quite rea- 
sonable given the well known insulating nature of alkanes 
and alkanethiols. 

While the solid curves in Fig. [3] are for a separation 
of 4.33A between the STM tip atom and the gold film, 
the (black) dotted curve for C6H12S is for a 5.05A sep- 
aration at which the calculated transmission Tfilm be- 
tween the STM tip and gold film is lower by a factor of 
approximately 35. Notice that the calculated values of 
/BEEM//tip are almost the same for the two separations, 
illustrating the insensitivity of the ratio /beem //tip to 
the separation between the STM tip and gold film in the 
tunneling regime, as has been discussed above. 

The calculated current ratios /beem //tip for some rep- 
resentative unsaturated molecules are shown in Fig. 3] 
Transmission resonances due to molecular states associ- 
ated with the unsaturated C-C bonds in these molecules 
give rise to dramatic increases in the BEEM current at 
bias voltages (indicated by arrows) at which the Fermi 
level of the STM tip crosses the resonant molecular en- 
ergy levels. As a result, even for CioHigS (depicted in 
the upper inset of Fig. (2) the largest value of /beem //tip 
is predicted to be within an order of magnitude of the 
NBT. For CeHioS (similar to CiqHisS but with a shorter 
C chain and only one H atom bonded to C atoms 3 and 4) 
and for CgHn (ethylmethylbenzene depicted in Fig. [T|) 
/beem //tip is predicted to exceed the NBT for bias volt- 
ages above the threshold for resonant transport4^ Thus 
experimental observation of single buried molecules such 
as these in BEEM experiments should be feasible with 
conventional BEEM equipment. Because the predicted 
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current ratios /BEEivi/J'tip for saturated molecules (Fig. 
[3]) are much smaller than those for unsaturated molecules 
of similar length (Fig. (H) it should be possible to study 
individual buried unsaturated molecules by BEEM and 
BEES by including a few unsaturated molecules in a 
self-assembled monolayer of saturated molecules of sim- 
ilar size on silicon under a thin gold film, the satu- 
rated molecules behaving as an insulating background 
surrounding the unsaturated conducting molecules. I.e., 
BEEM and BEES experiments extending to buried sin- 
gle molecules the classic STM studiesii^ of the proper- 
ties of single unsaturated molecules embedded in self- 
assembled monolayers of saturated molecules on gold 
substrates should be feasible. However, in interpreting 
such experiments it is important to note that the cal- 
culated current ratios /BEEM/-ftip shown in Fig. [3] for 
unsaturated molecules rise with increasing bias from val- 
ues far below the NBT. Because of this, experimentally 
observed values of the threshold voltage for the onset of 
the BEEM current are likely to depend on the sensitiv- 
ity of the BEEM equipment unless equipment capable of 
detecting extremely weak BEEM currents is used. 

As in Fig. [31 the solid curves in Fig. [Jare for a separa- 
tion of 4.33A between the STM tip atom and the gold film 
while the (black) dotted curve for CeHioS is for a 5.05A 
separation at which the calculated transmission Tfilm is 
again lower by a factor of approximately 35. Once again 
the calculated values of /BEEM/-^tip are almost the same 
for the two separations, illustrating the insensitivity of 
the ratio /BEEM/-^tip to the separation between the STM 
tip and gold film in the tunneling regime, this time for a 
resonantly transmitting molecule. 

In the course of the present work, preliminary results 
have also been obtained for some defects that may oc- 
cur in molecular SAMS on silicon covered with gold: It 
was found that small interstitial gold clusters embedded 
between decanethiolate molecules of the SAM and also 
single gold atoms substituting for H atoms of decanethi- 
olate molecules should give rise to BEEM transmission 
resonances and associated features in /BEEM/-^tip qual- 
itatively similar to shown in Fig. 31 However, these 
resonant features in Tbeem and Jbeem were found to 
be much weaker for the defects that were studied (a 6 
atom Au cluster roughly equidistant from the Au film 
and Si substrate and a single gold atom substituting for 
a H atom on the 5th C atom [from the Si substrate] of 
decanethiolate) than those of the unsaturated CioHigS 
molecule. 

Finally, it is interesting that the recent experimen- 
tal studiesi^ of gold-covered octanedithiol molecules on 
GaAs substrates observed BEEM current-voltage char- 
acteristics qualitatively similar to the resonant behavior 
in Fig. 31 However, since the experiments were on satu- 
rated molecules and the BEEM signal was observed only 
in a few very small patches of the sample the origin of 
the observed behavior is unclear. 

IV. CONCLUSIONS 

A better understanding of the electronic and structural 
properties of buried molecular layers and single molecules 



is essential for the development of molecular nanoelec- 
tronics. In this paper the feasibility of imaging single 
buried molecules and measuring their electronic spec- 
tra has been explored theoretically. It was found that 
the insulating nature of single alkane and alkancthiolate 
molecules chemisorbed on silicon substrates and covered 
with gold should make them difficult to detect by bal- 
listic electron emission miscroscopy and spectroscopy for 
all but very short molecules of this kind. However BEEM 
currents resonantly transmitted through molecules con- 
taining double C-C bonds or aromatic rings should be 
large enough to be accessible to present day BEEM equip- 
ment, and predictions of the BEEM current- voltage char- 
acteristics for some examples of molecules of each type 
have been presented. The predicted BEEM current- 
voltage characteristics for unsaturated molecules exhibit 
observable features (marked by the arrows in Fig. |4]) 
due to resonantly transmitting electronic states of the 
extended molecule that lie above the Fermi level of the 
metal overlayer. The energy differences between these 
resonant electronic states and the Fermi level of the metal 
film are given by \eV\ where V is the bias voltage be- 
tween the STM tip and metal film at which such a fea- 
ture occurs. Thus BEEM current-voltage characteristics 
can be used to measure the electronic spectra of buried 
unsaturated molecules in the environment between the 
metal and semiconductor layers. Since no bias voltage is 
applied between the metal overlayer and semiconductor 
substrate that surround the molecule these BEEM spec- 
tra differ from those obtained using conventional diode 
measurements in which the bias voltage is applied across 
the molecule and therefore modifies the electronic struc- 
ture that is being measured. 

At the present time rigorous theoretical techniques for 
molecular electronics, and especially molecular transport 
calculations, have yet to be developed. Thus the present 
study relies on a combination of semi-empirical and ab- 
initio methods as well as intuitively reasonable approxi- 
mations and assumptions. The results presented should 
facilitate the design and interpretation of experiments by 
providing theoretical answers to some rather basic ques- 
tions regarding BEEM and BEES of molecules that have 
not been addressed previously. Further theoretical work 
on this topic as well as experiments testing the present 
predictions are clearly desirable, and it is hoped that the 
present exploration will stimulate such efforts. 
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FIG. 1: (Color on line). The extended molecule studied in the BEEM transport calculations for the case where the molecule 
is ethylmethylbenzene." The STM tip is represented by 10 Au atoms in a (111) geometry and terminates in a single atom. 
The Au film is represented by a cylindrical cluster of 591 Au atoms in a bulk crystal geometry that presents (lll)-oriented 
facets to the STM tip and the molecule. The substrate is represented by an approximately hemispherical crystallite of 390 Si 
atoms. The Si dangling bonds are passivated with H atoms and a flat (111) surface faces the molecule. The molecule bonds to 
the Si substrate through a single C-Si covalent bond. A large number of single channel ideal leads (most of them not shown) 
representing the source and drain electrodes are attached to the atomic valence orbitals of appropriate surface atoms of the 
clusters representing the STM tip, Au film, silicon substrate and passivating H atoms. 



FIG. 2: (Color on line). Calculated electron transmission probabilities Tbeem{E , V) from electron source to drain via STM tip, 
gold film, molecule and and silicon substrate vs. electron energy E. Solid (dotted) curves are for at zero (2V) bias V between 
STM tip and gold film. The gold film Fermi energy is at i5 = and the Si conduction band edge at the Si surface Ec is shown 
by the dashed vertical line. The lower plots are for the decane thiolate C10H20S molecule in the lower inset. The upper plots 
are for the CioHigS molecule in the upper inset that exhibits resonant transmission due to the presence of a double C— C bond. 



FIG. 3: (Color on line). Calculated BEEM to film current ratios vs. tip bias voltage V for some saturated molecules: Decane 
(C10H21), decanethiolate (C10H20S) and hexanethiolate (C6H12S). The red horizontal dashed line is the nominal BEEM 
sensitivity threshold. Solid curves are for a separation of 4.33A between the STM tip atom and the gold film. The (black) 
dotted curve for C6H12S is for a 5.05A separation. 



FIG. 4: (Color on line). Calculated BEEM to film current ratios vs. tip bias voltage V for some unsaturated chain molecules 
CioHigS (shown in the upper inset of Fig. [2| and CeHioS, and for ethylmethyl benzene CgHn shown in Fig. [T] Arrows 
indicate values of the bias voltage at which the STM tip Fermi level crosses the energies of molecular resonant states. The red 
horizontal dashed line is the nominal BEEM sensitivity threshold. Solid curves are for a separation of 4.33A between the STM 
tip atom and the gold film. The (black) dotted curve for CeHioS is for a 5.05A separation. 
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